An interacting model for the cosmological dark sector 
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We discuss a new interacting model for the cosmological dark sector in which the attenuated 
dilution of cold dark matter scales as a"' i f{a), where f(a) is an arbitrary function of the cosmic 
scale factor a. From thermodynamic arguments, we show that /(a) is proportional to entropy source 
of the particle creation process. In order to investigate the cosmological consequences of this kind of 
interacting models, we expand f(a) in a power series and viable cosmological solutions are obtained. 
Finally, we use current observational data to place constraints on the interacting function /(a). 
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I. INTRODUCTION 

A better understanding of the physical mechanism be- 
hind the current cosmic acceleration is one of the ma- 
jor challenges both for the cosmology and fundamental 
physics. In the framework of the general theory of rela- 
tivity, the cosmic acceleration can be explained either if a 
new hypothetical energy component with negative pres- 
sure (dark energy) dominates the current composition of 
the cosmos [l[ or if the matter content of the universe is 
subject to dissipative processes (see, e.g. [2|). 

Following the first route, the simplest and most natural 
candidate for dark energy is the cosmological constant, 
A, which corresponds to the energy density stored in the 
true vacuum state of all existing fields in the Universe. 
From the observational point of view, it is well known 
that the ACDM model is in good agreement with almost 
all sets of cosmological observations. Despite its obser- 
vational successes, it suffers at least from two problems. 
First, and possibly the most serious one is the cosmolog- 
ical constant problem (CCP) in which the cosmological 
upper bound (pa < 1CT 47 GeV 4 ) differs from theoretical 
expectations (pa ~ 10 n GeV ) by more than 100 orders 
of magnitude. The other is known as coincidence prob- 
lem and consists to understand why pa is not only small, 
but also of the same order of magnitude of the energy 
density of cold dark matter (CDM) Q. 

An attempt to alleviate the cosmological constant 
problems is to allow the vacuum energy and dark matter 
to interact, leading to the so-called vacuum decay mod- 
els for which A is necessarily a time-dependent quan- 
tity. Cosmological scenarios with a dynamical A term 
were independently proposed about two decades ago in 
Refs. [ij (see also [H-IH for interacting models in the which 
the dark energy is represented by a smooth component 
parametrized by an equation of state pde = wpuE with 



w < 0). Most of these scenarios have specific decay laws 
for A which lead to particular cosmological solutions and 
phenomenology (see, e.g., Table I of @ for some proposed 
A decay laws). 

In this paper, we discuss a general cosmological sce- 
nario of A-CDM interaction in which the attenuated di- 
lution of CDM is described by an arbitrary function of the 
cosmic scale factor f(a). We show that this class of mod- 
els has many of the previous phenomenological scenarios 
as a particular case. Using thermodynamic considera- 
tions, it is shown that the function /(a) is proportional 
to the entropy source of the associated particle creation 
process. We discuss some interesting cosmological con- 
sequences of this A-CDM interaction, such as a possible 
loitering phase in a spatially flat universe and carry out a 
joint statistical analysis with recent observations of type 
la supernovae, baryonic acoustic oscillations and cosmic 
microwave background data to test its observational vi- 
ability. We work in units where (8ttG) 1/ ' 2 = c = 1. 



II. A-DARK MATTER INTERACTION 

According to the Bianchi identities, the Einstein field 
equations, G^ u = T^ v + Ag^, implies that A is neces- 
sarily a constant either if the energy-momentum tensor 
of matter fields and CDM particles T^ v is null or if it is 



separately conserved, i.e., UnT^ 



0. This amounts to 



saying that the presence of a time-varying cosmological 
term results in a coupling between A and other cosmic 
component. By considering the Friedmann-Lemaitre- 
Robertson- Walker space-time and a coupling between A 
and CDM particles, we have that 



Pdm + 3 — Pdm — ~PA — —Q , 

a 
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where pdm and p\ are the energy densities of CDM and 
A respectively and the dot sign denotes derivative with 
respect to the time. 

As the cosmological term interacts with CDM parti- 
cles, their energy density will dilute slowly compared with 
the standard relation, pd m oc a~ 3 . Thus, the deviation 



2 



from the standard dilution may be characterized by the 
function /(a), such that 



Pdm 



By substituting the relation @ into (P), we find 



dpA 
da 



da ' 



(2) 



(3) 



where for f(a) = const. = pd m ,o the standard ACDM 
scenario is fully recovered. Note that for f(a) = pdm,o^ e 
the dynamical A(i)CDM scenarios recently discussed in 
Refs. [IM1 are obtained whereas the vacuum decay 
model with A cc H investigated in Refs. [l4| is recov- 
ered when f(a) = c\ + C2<i 3 / 2 . More general coupled 
quintessence scenarios in which the coupling is written 
in terms of the scalar field <f> as Q = y^2/ 3(3 pd m (f> lh\ 
(see also [16|), where (3 is a coupling function, are recov- 
ered by identifying 



/ 



(4) 



In this paper, we will work with the /(a) function as 
a power series of scale factor, i.e., 



f(a) = fo + fia + f 2 a 2 + 



(5) 



where fo, fi and f 2 are constants. From Eqs. ©, §5$ 
and (0, we find 



Pdn 



f a- 3 + ha- 2 + ha- 1 + 



h _ 2 

Pa = Pa,o + y a + 2/ 2 a 



(6) 



(7) 



where pa,o is & n integration constant. 

Now, neglecting the radiation contribution and consid- 
ering that the baryonic content is separately conserved, 
the expansion rate of the Universe can be written as 



3 

n = nA,o + (tt +tt b . )a- 3 + -n 1 a- 2 - 



(8) 



where H = (H/H ) 2 , f2 A ,o = PA,o/p c ,o, ^fc : o = Pb : o/p C: o 
represents the density parameter of baryons, f^o = 
fo/p c ,o, = h/pc,o (with i = 1, 2, n) and p cfi = 
is the present critical energy density. In what follows, we 
will restrict our analysis up to the first order expansion 
of f(a). The practical reason for this choice is shown 
in Figure 1 where the ratio pa/ Pdm as a function of the 
redshift is displayed for first (solid line) and second order 
expansions in /(a) assuming the best-fit values of Qq, 
fli and f2 2 discussed in Sec. IV. Clearly, higher order 
terms do not modify the ratio between A and the CDM 
component for the interval considered. We note that an 
expansion rate similar to Eq. ([SJ has been explored in 
reconstruction studies of the dark energy density (see, 
e.g., Ref. [13). 




FIG. 1: The ratio of the energy densities as a function of the 
redshift. It was used the best-fit values of f2o, ^i and fl2- 



III. THERMODYNAMICS OF INTERACTION 

Let us now consider some thermodynamic aspects of 
interaction between the vacuum and CDM particles. The 
thermodynamic behavior of a decaying vacuum system is 
simplified either if one assumes that the chemical poten- 
tial of the vacuum component is zero or if the vacuum 
medium plays the role of a condensate carrying no en- 
tropy, as happens in the two fluid description employed 
in the super fluid thermodynamics [l8j . 

In this case, the thermodynamic description requires 
only the knowledge of the particle flux, N a — nu a , and 
entropy flux, S a = nau a , where n = N/a 3 and a = 
S/N are, respectively, the concentration and the specific 
entropy (per particle) of the created component. The 
equation for n has necessarily a source term, i.e., 



3— n 

a 



* = nT 



(9) 



where 'J is particle source ("J > 0), or a sink (vt < 0), 
and r is the rate of change of particle number. Since 
Pdm, = nm, Eq. © gives 

n = n Q a- 3 f(a) . (10) 

Substituting above equation in Eq. ©, it follows that 



/ = /r. 



(ii) 



As argued in Ref. [18| , the energy exchange between the 
vacuum and the cold dark matter may take place in sev- 
eral ways. The most physically relevant case has been 
termed adiabatic decaying vacuum. In this case, several 
equilibrium relations are preserved, and perhaps more 
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which includes 557 data points after selection cuts. 
Following Ref. [22| we also use measurements derived 
from the product of the CMB acoustic scale i A = 
i:d A {z if )/r s {z Sf ) and from the ratio of the sound hori- 
zon scale at the drag epoch to the BAO dilation scale, 
r s (zd)/ Dv(zbao), where d A (z*) is the comoving angular- 
diameter distance to recombination (z* = 1089) and 
r s (z*) is the comoving sound horizon at photon decou- 
pling. In the above expressions, Zd — 1020 is the red- 
shift of the drag epoch (at which the acoustic oscilla- 
tions are frozen in) and the dilation scale, Dy, is given 
by Dv(z) = [zr 2 (z)/H(z)] 1 ^ 3 . By combining the ratio 
r s {z d = 1020)/r s (z* = 1090) = 1.044 ± 0.019 013 
with the measurements of r s (zd)/ Dv(zbao) a t zbao — 
0.20, 0.35 and O.6., one finds / . 20 = d A (z*)/D v (0.2) = 
18.32 ±0.59, /0.35 = d A (z*)/D v (0.35) = 10.55 ± 0.35 , 
fo.e = d A (z*)/D v (0.60) = 6.65 ± 0.32 JHHi]. 

In our analysis, we minimize the function \\ — XsNe + 
Xcmb/bacp where xInc and Xcmb/bao correspond to the 
SNe la and CMB/BAO \ 2 functions, respectively. Fig- 
ure 2 shows confidence regions (68.3% CL and 95.4%) 
in SIq — ill plane obtained from this joint analysis. At 



FIG. 2: The results of our statistical analysis. Contours of x 2 2er we found f2 = 0.25±qq3 and fii 



-0.028 
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in the plane Qo 
2.30 and 6.17. 



fli. These contours are drawn for Ax 



importantly the entropy of created particles increases 
wheras the specific entropy (per particle) remains con- 
stant [a = 0). This mean that 



S 
S 



N 
N 



(12) 



From Eqs. (jlll) and (TT2"j) it is straightforward to show 
that 



f<xS. 



(13) 



The second law of thermodynamics, therefore, implies 
that /(a) must be an increasing function of the scale fac- 
tor, which requires fli > 0. Physically, this result indi- 
cates that only an energy flow from vacuum to cold dark 
matter is allowed. As discussed in the next section, this 
is in disagreement with current observational constraints 
that favor an energy flow from dark matter to vacuum 
(fil < 0). 



IV. OBSERVATIONAL ANALYSIS 

In order to delimit the parametric space ilo — 0,% we 
perform a joint analysis involving current SNe la, BAO 
and CMB data. In our analysis, we fix Jl^o = 0.0416 
from WMAP results [lj| (which is also in good agreement 
with the bounds on the baryonic component derived from 
primordial nucleosynthesis [2fjj ) and marginalize over the 
Hubble parameter Hq. 

We use one of the most recent SNe la compilation, 
the so-called Union 2 sample compiled in Ref. [2lJ 



expected, we found values of fij ~ since it acts as a cur- 
vature term (a~ 2 ) in Eq. ©. As shown in Section III, 
the interacting parameter f^i is restricted from thermo- 
dynamical arguments to be positive. By considering this 
physical constraint on Qi and minimizing the function 
Xt, we obtain Q = 0.24±°;°| and fii = 0.0±g$ at 2a 
level. An observational comparison between the model 
here described and the coupled quintessence scenario of 
Ref. [l5[ can be made by comparing the above estimates 
for fii with the constraints on the coupling function /3 
[see Eq. Q]. For the case in which j3 remains con- 
stant with z, Ref. [26| used LSS and CMB data to found 
|/3 1 < 0.1, which is in full agreement with the bounds on 
f2i obtained in our analysis. 

For the sake of completeness, we also derive the time 
evolution of the density parameters Clb(a), £lo(a), fli(a) 
and Qa(o>) which are given, respectively, as 

f\ a~ 3 n , x 1.5ftiar 2 f2 A ,o 
S2j(a) = —. . , ih(a) = — — — — , $2 A (a) 



C(«) 



C(«) 



where 



£(a) = fi A , + (O + fi 6 , )fl 3 + 2^i«" 2 

and i = 6,0. Figure 3 shows the evolution of Q,h(a), 
flo(a), Oi(a) and ^Ia(o) with the logarithm of the scale 
factor ln(a) for f7 h , = 0.0416 and fl = 0.22, 0.25, 0.28 
and fii = -0.078, -0.028, 0.022, which correspond to 
lower limit, best- fit and upper limit (at 2a of confidence) 
from statistical analysis described above. In all cases 
a mix of baryons (< 20%) and dark matter (> 80%) 
dominates the past evolution of the Universe whereas the 
contributions relative to the integration constant Oa.o 
and the interaction parameter fii vanish at high-z. As 
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FIG. 3: Evolution of the density parameters as a function of ln(a) for the lower limit (left panel), best-fit (central panel) and 
upper limit (right panel) on Qq and Qi obtained from the statistical analysis described above. 



expected, and in agreement with observational data, we 
note that £l\,o is the dominant component from a value 
of a* < 1 on, while the interaction parameter Qi(a) is 
always the subdominant throughout cosmic evolution. 

Another interesting aspect of this fiat vacuum decay 
model is that it can accommodate a loitering phase. As 
discussed in Ref. |27|, a loitering phase is characterized 
by a dH/dz = 0. Deriving Eq. (JSJ and using the loitering 
condition, we hnd (1 + z = a -1 ) 



In the context of this flat interacting model, the ex- 
istence of a recent loitering phase is possible only if 
|fii| > CIq + £lb,o- By considering the obsevational limits 
on these parameters given above, we note that a realistic 
loitering phase is possible only in the future (zi < 0). 
Finally, by considering that a /a = H + H 2 , it follows 
that (a/a) z=Zl = H 2 (z = zi), which clearly shows that 
models with a loitering phase are compatible with an 
accelerating universe, as indicated by SNe la data. 



/(a) is proportional to the entropy source of the particle 
creation process, which means that only an energy flow 
from vacuum to cold dark matter is allowed. We have 
investigated the observational viability of this class of 
models from recent data of SNe la and the distance ratio 
from baryon acoustic oscillation at different rcdshifts and 
CMB decoupling at zls = 1089. In particular, when /(a) 
can be expanded in a power series of the scale factor, we 
find strong constraints on the parameters of this model, 
with a slight deviation from standard ACDM dynamics. 

Finally, it is worth emphasizing that we have restricted 
the present analysis to coupled quintessence models in 
which ?«0 = —1 (dynamical A models, whereas a full 
treatment of the dark matter-dark energy interaction 
must also take into account the role of the dark energy 
equation-of-state in the process. Some theoretical and 
observational consequences of a w-CDM interacting sce- 
nario with a time-dependent coupling term, as well as a 
scalar field description for this class of models is currently 
under investigation. 



V. CONCLUSIONS 

We have discussed cosmological consequences of a new 
cosmological scenario of vacuum decay in which the de- 
viation of standard evolution of the cold dark matter 
component is characterized by an arbitrary function of 
the cosmic scale factor /(a). Using thermodynamic ar- 
guments, we have shown that the interaction function 
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